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Abstract 

We study the properties of neutron stars adopting relativistic equa- 
tions of state of neutron star matter, calculated in the framework of 
the relativistic Brueckner-Hartree-Fock approximation for electrically 
charge neutral neutron star matter in beta-equilibrium. For higher 
densities more baryons (hyperons etc.) are included by means of the 
relativistic Hartree- or Hartree-Fock approximation. The special fea- 
tures of the different approximations and compositions are discussed 
in detail. Besides standard neutron star properties special emphasis 
is put on the limiting periods of neutron stars, for which the Kepler 
criterion and gravitation-reaction instabilities are considered. Fur- 
thermore the cooling behaviour of neutron stars is investigated, too. 
For comparison we give also the outcome for some nonrelativistic equa- 
tions of state. 
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I Introduction 



A necessary ingredient for solving the structure equations for (rotating) neu- 
tron stars (NS) is the equation of state (EOS) [Q. For NSs the EOS has to 
cover a wide range of densities ranging from super-nuclear densities (up to 
5 to 10 times normal nuclear matter density) in the star's core down to the 
density of iron at the star's surface. At present, neither heavy-ion reactions 
nor NS data are capable to determine the EOS accurately, and the behaviour 
of high-density matter is still an open and one of the most challenging prob- 
lems in modern physics, containing many ingredients from different branches 
of physics. The theoretical determination of the EOS over such an enormous 
range has therefore to rely mainly on theoretical arguments and extrapola- 
tions for which no direct experimental confirmation exists. The best one can 
do in such a situation is a step-by-step improvement of the available models 
for the EOS. Since the central density of a NS is so extreme, both the Fermi 
momenta and the effective nucleon mass are of the order of 500 MeV, one 
should prefer a relativistic description 0-0. 

Neutron star matter differs from the high density systems produced in 
heavy ion collisions by two essential features: a) Matter in high energy col- 
lisions is still governed by the charge symmetric nuclear force while neutron 
star matter (NSM) is bound by gravity. Since the repulsive Coulomb force is 
much stronger than the gravitational attraction, NSM is much more asym- 
metric than "standard" matter, b) The second essential difference is caused 
by the weak interaction time scale of ~ 10~^'^s, which is small in comparison 
with the lifetime of the star, but large in comparison with the characteris- 
tic time scale of heavy ion reactions. For that reasons "normal" matter is 
subject to the constraints of isospin symmetry and strangeness conservation, 
but NSM has to obey the constraints of charge neutrality and generalized 
beta-equilibrium with no strangeness conservation. It is obvious from these 
considerations that NSM is an even more theoretical object than nuclear 
"normal" matter, with a rather complex structure 

Due to these features one can adopt the following structure of a neutron 
star: The two outer crusts of the star have crystalline structures, for which 
rather reliable EOSs exist in the literature. In the uniform density region 
(p > xlO^^ g/cm^) of neutron star matter, one has to deal with a system of 
interacting baryons (i.e., p, n. A, S, possible A's, etc.) and/or quarks {u,d,s 
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flavours, uniform or not uniform) that are in generalized /3-equilibrium with 
leptons (e~,/i~) [Q-JITI. Furthermore one may encounter meson condensates 

[0, 0, 0- 

Despite the differences of hot symmetric non-strange matter produced in 
high energy collisions and cold asymmetric charge neutral and strangeness 
containing matter of NSs one should combine both systems in a common the- 
ory. This is possible to a certain degree in modern field theoretical relativistic 
approaches, and NSs are unique systems in the sense, that they offer a test 
bench for the EOS of exotic NSM, which can not be mimicked in terrestrial 
laboratories. 

The first attempt, based on the Fermi gas model, to incorporate the role 
of the hyperons in NSs is due to Ambartsumyan and Saakyan ||13|, which 
was later improved in the nonrelativistic scheme by several authors [1^. The 
first systematic investigation in the relativistic framework was performed by 
Glendenning, who used the relativistic mean-field approximation with in- 
clusion of hyperons and A's. 0. This model was later extended by using 
the relativistic Hartree-Fock-approximation (RHF) [|^, |^, ^ and by improv- 
ing the Lagrangian in the mean-field approximation (or relativistic Hartree- 
approximation; RH) |Tl], |T5|, |T6|. One of the most exhaustive investigations 
in the later framework was performed by Schaffner and Mishustin, who used 
modern phenomenological interactions in the nucleon and hyperon sector, 
respectively. In this last model kaon condensation turned out to be unlikely 

From a microscopic standpoint such treatments are not satisfactory, since 
the interaction is adjusted in a phenomenological manner to properties of fi- 
nite nuclei and the parameters of the Bethe-Weizsacker mass formula. Other 
drawbacks are, for instance, the vanishing of the vr-meson contributions in the 
mean-field approximation, and large a-meson self-interactions are needed to 
reduce the incompressibility K. Some of these deficiencies do not occur in 
the framework of the relativistic Hartree-Fock approximation, which also 
resembles, in both its mathematical structure and its Lagrangian density 
rather closely to the microscopic relativistic Brueckner-Hartree-Fock the- 
ory, for which one-boson-exchange potentials (OBEP) are adjusted to the 
two-nucleon data. In order to incorporate a more microscopic description of 
neutron star matter, we have recently developed a model in which this mat- 
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ter is described for moderate densities by self-consistent RBHF-calculations 
for /3-stable matter (n, p, e~ and fj,~ only) [0, [T^ . Due to technical diffi- 



culties, RBHF-calculations are presently restricted to densities up to 2-3 
times nuclear matter density. Since at higher densities other baryon states 
become populated and RBHF-calculations with inclusion of such states are 
presently not feasible, we extrapolated the EOS at higher densities within 
the RH and/or the RHF approximation. The essential new feature of this 
scheme in comparison with older treatments was the different adjustment of 
the parametrization of the RH and RHF Lagrangians. Normally one restricts 
oneself to the reproduction of the properties of symmetric nuclear matter, 
but in the new parametrization we used the outcome of RBHF calculations 
of asymmetric and NS-matter as well. It turned out that the RH approxi- 
mation is not so flexible in reproducing the properties of NSM (asymmetry, 
composition etc.). Since the Hartree approximation has been used in the 
vast majority of the earlier investigations, we will include this approximation 
in our considerations, too, which has also the advantage of greater trans- 
parency. Special attention will be paid to the RHF-approximation, where 
new features enter into the properties of NSM. In both approximations it is 
possible to incorporate more baryon states (for more details, see Refs. @, 0, 
and following sections). 

In order to achieve a certain degree of self containment we recapitulate also 
the basic theory for the EOSs, and for the properties of neutron stars, i.e. 
structure of rotating stars in general relativity, stability of rotations (Kepler 
criterion, gravitation-radiation instabilities), cooling of neutron stars etc. 
The contribution is organized as follows. Section II is devoted to the EOS of 
NSM, where we discuss, after a brief theoretical review the different models 
of the EOSs. The next section deals with the properties of NSs and the 
summary is given in Section IV. 



II Equation of state 

a) General theory: 



The EOSs of neutron star matter were determined as described in Refs. [|T7 
18| . For the two outer crusts we used EOSs taken from the literature [0, ^ 



The general dynamics of the hadron/lepton system in the uniform region of 
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neutron star matter is governed by an OBE Lagrangian of the following form: 
C{x) = J2 ^b(^) (II.l) 

B=p,n,S±o,A,HO.-,A-.o.+>++ 




B=p,n,...A 




The forces between the different baryons are mediated by the exchange of 
different mesons M = a,uj, . . .The leptons e~ and /i~ are treated as free 
particles. Furthermore one has to impose the constraints of beta-equilibrium 
[qb and fiB denote the electric charge and chemical potential of the baryon 
B, respectively): 

fJ'B = f^n- qBfJ-e , fJ'p = fJ'e ■ (11.2) 

and charge neutrality. For the many-body treatment we employed the rel- 
ativistic many-body Green's function scheme. Here one has to solve - for 
Brueckner-type approximations - a coupled system of equations which con- 
sists of the Dyson equation for the Green's function G, the effective scattering 
matrix T in matter, and the equation for the self-energy S of a baryon in 
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(G°)-i(l,2)S(l,2)G(2,l') = 5(1,1') (II.3) 
<12|r|l'2'> = < 12|V|1'2' - 2'1' > + (II.4) 
i < 12|V^|34 > A(34,56) < 56|T|1'2' > , 

with 

S(l, 2) = - < 14|T|52 > ^(5, 4) . (II.5) 

For the intermediate baryon-baryon propagator A, we have chosen in 
the RBHF approximation the Brueckner propagator. V denotes the OBE 
potential, which has the following structure 

< 12|V|1'2' > = XI < 12|Va/|1'2' > . (II.6) 

M=cr,w,p... 

The RH and RHF approximations are obtained for T = V and T = V — V^^ , 
respectively. For the OBE potentials we selected the potentials A and B 
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constructed by Brockmann and Machleidt ||2^. In the RBHF-treatment, 
performed in the full Dirac space, of symmetric and asymmetric nuclear 
matter both potentials (they differ in the strength of the tensor force, which 
increases from A to B) give good agreement with the nuclear matter pa- 
rameters {E / A, pQQ, K^, J), and also the volume parameters L and Kgym of 
asymmetric matter are in accordance with the data. An illustrative compar- 
ison with other treatments is given in Table I ||2ll . 



Up to 2-3 times equilibrium nuclear matter density we used for the EOSs 
the outcome of an RBHF-calculation for neutron star matter by restriction 
to p,n,e~, and only. For higher densities more baryon states B become 
populated in /3-stable neutron star matter. For high density systems with 
such a complex composition the RBHF scheme is not feasible. In order to find 
an extension of the microscopic RBHF-EOS to higher densities we selected 
both the relativistic Hartree and the relativistic Hartree-Fock approximation, 
where the parameters in the nucleonic sector (coupling constants etc.) are 
adjusted to the results of RBHF-calculations for asymmetric and NS matter. 

Explicitly, the Lagrangian density for the RHF-approximation, where 
the forces are mediated by the exchange of a-, u-, and p-mesons, and a 
self-interactions are included, is given by |T7[ : 



B 



Anir 



I'D 



fix CT^^" 
— I'l'TB ■ d,7T - gf^T ■ p^{x) - ff- T ■ F 



+ i [d^aix)d''cr{x) - mla^ix)] + ^ [d^7z{x) ■ d''7z{x) - mliT^x 



- -^F';,{x) ■ F^^'^ix) + ^myix) ■ p^ix) - ^F-,ix)F^^'-ix) + ^m>^(a;)^,(x 
1 . . , 1 



^mNbN[ga(^ix)f - -CN[gaCrix)]^ , (11.7) 



with 



F^A^) = df,uj^{x) - d^ujf,{x) , F''^^{x) = d^p^{x) - d^p^{x) . (II.8) 

With respect to the saturation properties of infinite nuclear matter (INM) 
(see Table I) the results are identical to the RBHF-outcome. It seems that, at 
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present, such a procedure is the only possibihty to incorporate more baryons 
and to estabhsh a connection to microscopic RBHF-calculations with re- 
ahstic OBE-potentials simultaneously. In this context one has to remark 
that the standard relativistic approach in almost all cases is based on a 
pure phenomenological RH-treatment, where the parameters in the nucle- 
onic sector are adjusted to the properties of symmetric nuclear matter, so 
that our approach has the advantage to contain more microscopic elements 
(for more details, see Refs. [l^, ^). In the calculations we used an im- 



proved parametrization [|T8|, and the parameter sets are given in Table II. 



b) Comparison of the different approximations 

First one has to test whether the described approximations can reproduce 
the properties of asymmetric and NS matter in the density range, where the 
RBHF-treatment is applicable (nucleons and leptons only). In Fig. 1 we 
compare the EOSs for different asymmetries in the different approximations. 
The agreement of the RBHF-EOSs with the RHF-EOSs is rather satisfactory 
for the whole asymmetry range. However, for the RH-EOSs the agreement 
is worse for larger asymmetries, furthermore the RH-EOS becomes stiffer 
for higher densities. This behaviour of the RH-approximation can better 
be inferred from Fig. 2, where the comparison for the pressures for NSM is 
shown. In Fig. 3 we exhibit, as a further example, the comparison with re- 



spect to the baryon/lepton composition. As discussed in Refs. [|T^, |T8| the 
RHF-approximation with no p-tensor coupling gives the best agreement with 
the RBHF-EOS. The RHF-approximation has the additional advantage that 
the Lagrangian density and the mathematical structure resembles more the 
RBHF-scheme than the RH-approximation (no 7r-meson, exchange contri- 
butions etc.). For the presentation we selected the Brockmann-Machleidt 
potential B [^]; for the potential A the situation is completely analogous 

mm- 

c) Neutron star matter 

a) General considerations 

For the calculation of NS-properties the EOS is needed for a wide range 
of densities, stretching to several times of nuclear matter saturation density 
[Q]-0. If one extrapolates the described scheme to the density domain of 
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NSM, one faces the following dilemma: It is well known that NS-properties 
depend strongly on the properties of the EOS near saturation, which is obvi- 
ous for lighter NSs but also true for heavier stars p6[. Therefore one should 



use in this region an EOS, which is either based on microscopic RBHF- 
calculations |2l|, ^ or phenomenological parametrizations of the 

RH/RHF-approximation, adjusted to nuclear data [23|. In both cases one 



obtains inevitably a rather stiff EOS caused by the low value of the Dirac 



mass of about 0.6 m^r (correct spin-orbit splitting [29|; the RBHF-scheme 
gives also reasonable results for finite nuclei |Q) at saturation p3|, |3l|, ^ 



The resulting meson fields are then rather large and consequently one ob- 
tains a sharp drop of the Dirac masses with increasing density. This fea- 
ture is even amplified by the unavoidable occurence of negative values for 
CN{bj\f/c]y ~ —1) 0, which causes a nonmonotonic behaviour of the 

effective cr-mass 



*2 



ml + gl \bNmN < QaO- > +Cn < QaO- 



(11.9) 



which increases the attraction beyond < >°= —hNrri]^ /2cn- As long 
as the composition of NSM is restricted to n,p, e~, and /i~ only, the EOS is 
sufficiently stiff to reach the necessary central pressure of the star at moderate 
densities (see Section III. e). However if one includes more baryons in the 
NSM-composition the EOS becomes considerably softer and higher densities 
are needed to obtain sufficient central pressure. The scalar fields are therefore 
in this case rather large and hence negative Dirac masses for the nucleons 



occur in the calculation |Tl|, [T^. One has tried to overcome this problem 
in a phenomenological manner using so-called stabilized a-functional forms 
0, for which the dangerous negative curvature of m*^ is switched off. 
However for the familiar parametrizations PL-2 and PL-40 negative Dirac 
masses also still occur in NSM. Schaffner and Mishustin have circumvented 
this problem by the rather questionable ad hoc assumption of using always 
absolute values for the Dirac masses, so implicitly changing the stiffness 
of the EOS [|lT|. Another approach by which one gets now along with a 
reduced attraction, i.e. bN,CN > 0{cn ^ ^a?) reduces the repulsion by an 
additional cu-self-interaction [03, 1351 . In this manner one obtains softer EOS 



with smaller cr-fields, where negative Dirac masses do not occur. However 
the resulting EOS may be too soft and an additional repulsion among the 
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hyperons seems necessary Furthermore the asymptotic behaviour is 

changed from the standard behaviour proportional to to that of an ideal 
gas (oc Pb^). Unfortunately the NSM-results for this case were not applied to 
NS-properties in Ref. [0]. For the sake of comparison we performed therefore 
some NS-calculations and obtained - as expected - maximal masses around 
1.5 Mq. (For instance, for the parameter set TM 1 we obtained 1.561 Mq 
with a central energy density ec = 740 MeV/fm^ llMl-) ^'^^ might wonder 
why this problem was not discussed in the other investigations of NSM. 
In the vast majority one prefers an EOS which is based on the standard 
nuclear matter parameters but with a Dirac mass of approximately 0.78 m^v 



[y, |g, |T5|, [Ty, ^ . Only in this window |32[ one obtains positive values for 
cn (or b]\r ^ \cn\) and more moderate meson fields, so that the problem of 
negative nucleon Dirac masses is avoided. The reasons given for this choice 
rest on a reproduction of the effective mass (~ 0.83 itin) |H, |36|, however a 
closer inspection according to a more elaborate investigation by Celenza and 
Shakin shows that values of again 0.6 mjy are more appropriate for the Dirac 
mass [0 (see also Ref. @]). Furthermore one favours rather low values for 
the saturation density and higher values for the incompressibility in order 
to stiffen the EOS (see, for instance, Ref. 0). For instance, according to 
the Hugenholtz-vanHove theorem, (yf^^v varies by approximately 40% in the 
range of 1.3 fm"^ < p^ < 1.42 fm"^ 



Common to all these approaches is to impose a certain behaviour of the 
NSM-EOS in the NSM-domain, which is not known. Controlled is this be- 
haviour by an additional parameter (large Dirac mass; switching off parame- 
ter; coupling constant of the vector self-interaction etc. H, |, ^ ^ ^ |35|). 



In order to overcome this dilemma, namely to keep the connection to the 
EOS in the vicinity of nuclear saturation and the described unpleasant fea- 
tures of the NSM-EOS with additional baryons one is also forced in our 
approach to implement a working extrapolation hypothesis controlled by an 
additional parameter. Also RBHF-calculations in symmetric matter favour 



higher Dirac masses in this density region |38|. In order to keep the scheme as 



simple as possible we extrapolate as follows: We maintain the dynamics till 
the maximum of the effective cr-mass and extrapolate beyond this point via 
a Lagrangian of the same structure, but with new self-interaction couplings 
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and a modified cr-mass: 



ml - = ml + ^mlh\ f- - l) , (H-IO) 
i>N — ^ i>'^ = a— ; cat — = a . (H-H) 

Cat 

For this model the effective cr-mass and its derivative agree with the mi- 
croscopic model at the transition point (see Eq. |11.9| ). a = gives the 
original dynamics; a = results in a linear Walecka model, which turns 
out to be insufficient |T^ (negative Dirac masses). With a > one can 
now control the stiffness. One should remark in this context that the main 
effect of the extrapolated dynamics is a more moderate drop of the Dirac 
masses beyond four times nuclear matter saturation density in accordance 
with Ref . p8| (see Fig. 4). With respect to the pressure, the differences 



against the original dynamics are rather small (10%) in the lower part of 
the NSM-density domain and become very small in the high-density NSM- 
region, where the unchanged cu-meson repulsion dominates. We have also 
tested the pure linear or quadratic extrapolations, i.e. b'jy = biy/2, c'^ = or 
h'^ = 0,c'j^ = —Cn', m'^ = m^, but their results can be incorporated closely 



in the scheme selecting special a-values . 

P) Results and discussion 

Essentially three main features characterize the properties of NSM: a) 
the stiffness of the EOS b) the relative coupling constants of hyperons, c) 
the chosen many-body approximation. All three points are correlated, but 
nevertheless we will try to separate them to a certain degree in order to 
extract some insight into the structure of the problem. 

The influence of the first point is obvious, since the necessary central pres- 
sure for stable stars can be reached earlier for larger stiffnesses. Important 
is the second issue. The relative hyperon couplings are not well known, since 
the information from hypernuclei data (for more details, see, Refs. [p9| , ^) 



permits a wide bandwidth (for instance, 0.4 < XHa = QHa/gNu < 0.8 in the 
relativistic mean field approximation). Therefore a number of choices were 
made in the literature, reaching from universal coupling, which gives a first 
insight, to ratios motivated by the quark model 0, |, ^ ||, ||, |ll|, |15|, [TBI. A 
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suitable first choice in the latter case is to use the SU(6) symmetry for the 
vector couplings i.e., 



-^Qu^N = ^9u;A = ^9ujT. = 9ojS ] 9pN = 9pT. = 9p=. , 9pA = 0, (11.12) 

and to fix the cr-coupling according to the potential depth of the hyperon- 
particle in nuclear matter (~ —30 MeV). According to these uncertainties 
we will investigate the problem in a more systematic manner by making 
several choices for the cr-hyperon coupling and fix the cu-coupling by means 
of the hyperon potential depth |]15|, The selected choices for the EOSs 
are described in Table III. The hyperon couplings are given in Table IV. 

The general trend is that the EOSs become softer with decreasing cou- 
plings, because then the conversion of nucleons into hyperons is energetically 
favourable, due to the smaller repulsive force dominating in the high-density 
domain of NSM. A special feature of the Hartee approximation are large p- 
couplings, which are necessary for the adjustment of the symmetry energy 
1^. For that reason the charge-favoured (isospin unfavoured) A~ does not 



occur (or play a minor role; see discussion of the RHF-approximation) and 



therefore is usually neglected a priori in this approximation [T^, |T5|, p!q] . 
For the S~ the charge compensation still dominates and so the S~ occurs 
rather early. Increase of the p-coupling, for instance hke = ^gp^, would 



reverse the onset of S and S+ |T3. A further and more severe point is 



the selection of the many-body approximation. For the Hartree approxi- 
mation we obtain compositions and stellar properties which are familiar to 
investigations performed earlier in this framework. This approximation is 
characterized as a high density approximation with a relatively stiff EOS 
containing no exchange contributions. For the latter reason the coupling 
constants are larger than in the RHF-scheme. Two examples for the com- 
position (universal coupling and SU(6) scheme) are given in Figs. 5 and 6, 
from which one can infer the dependence on the hyperon couplings. For the 
universal coupling the attraction is more favoured, which reduces the Dirac 
masses. Therefore the onset of the hyperons occurs earlier and also A's are 
possible. This causes a lower pressure at lower densities, however for high 
densities the stronger repulsion causes then a stiffer EOS (see Fig. 7). For 
these reasons one should obtain smaller NS-masses for the universal coupling 
for smaller central densities than for the SU(6)-coupling. For higher densities 
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the situation is reversed (see Sec. IlI.e). 

Due to the additional degrees of freedom the RHF-EOS is softer than 
the corresponding RH-EOS in a large density domain. A special new feature 
of the RHF-approximation is that the onset of a particular baryon species 
in NSM is, as in the RH-approximation, solely controlled by the Hartree 
term. Due to the lower couplings the participation of the other baryons 
is more enhanced and in general the onset of the hyperons occurs earlier 
than in the RH-scheme, resulting in a quite different baryon/lepton compo- 
sition and a softer EOS. A further characteristic difference is caused by the 
smaller p-coupling constants in the RHF-theory (see Table II). Responsible 
for this decrease of the coupling constants Qp are - as in the case of realistic 
OBE-potentials in the RBHF-theory - the exchange contributions, which 
contribute to the symmetry energy. For that reason the charge-favoured but 
isospin unfavoured A~ plays now - as in former nonrelativistic many-body 



approximations with correlations [Tj] - an important role in the composition. 
One might further wonder that one has not included the tu-tensor coupling 
for the hyperons, since, for instance, according to the quark model f/}/gl^ 
becomes —1 for the A hyperon. For that reason we included this term in a 
test calculation but the impact on the EOS was rather small and can 
therefore be neglected in comparison with the other uncertainties. 

If one would use naively the same hyperon couplings as in the RH- 
approach one obtains compositions - shown for the SU(6) in Fig. 8 - where 
the A's occur rather early (for universal coupling the baryons occur in the 
order E~,A, A~ with a high contribution of A's at higher density [|18]). For 
that reason the EOSs are rather soft at moderate densities and should cause a 
weak increase of the NS-mass in the lower part of the central energy-density 
region resembling in this part the behaviour of NS-masses calculated with 



EOSs with rather low incompressibilities [0, ^ (see Fig. 17). The strong 



abundance of the A's in such EOSs is caused, as discussed before, by the 
early onset of the charge-favoured A~, which is not so strongly isospin hin- 
dered in the RHF-approximation. The onset of hyperons in this model may 
be overestimated, since it is controlled by their Hartree contribution solely. 
A special feature of the RHF-theory seems to be that the Fock contributions 
for the scalar part of the self-energy approximately cancel each other but for 



the vector part they amount to 50 % of the total value [ffH]. For that reason 
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the hyperon couplings taken from the RH-approach do not reproduce the 
potential depths of the hyperons in nuclear matter in the RHF-scheme (at- 
traction too strong) and hence favour the hyperons. If one now corrects this 
deficiency by adjusting the hyperon couplings {x^^h (RHF) < x^^h (RH); see 
Table IV), the composition becomes now dominated by the A's (see Fig. 9). 
The occurence of the hyperons (A-hyperon) is rather late even for weak cou- 
plings (for instance, for Xujh = 0.5 at p ~ 0.5 fm~^). The peculiar behaviour 
of the pressure disappears now and consequently the NS-masses as func- 
tion of the central energy-density give now the standard strong increase for 
moderate densities (see Fig. 19, 20). The models treated so far involve the 
assumption that the A-coupling agrees with the nucleon coupling. But ac- 



cording to investigations of ter Haar and Malfliet [43|, the A-mass does not 



change very much from the vacuum case, indicating smaller A-couplings, 
which would delay the onset of the As. In a recent treatment by R. Rapp et 



al. the choice g^A = 9 an = 0.625 qnn was recommended, but one could 
use as an option also this choice for the attraction only [jl4|. For the first case 
one gets a slight shift of the A-threshold towards higher densities (see Fig. 10) 
and the EOS becomes much softer due to the lower repulsion. This can be in- 
ferred from Fig. 20, where the maximum star mass is approximately 1.5 Mq. 
In the last case one deals with a strong repulsion combined with moderate 
masses for the A's. For that reason one expects a minor role to be played by 
As. The calculation of the chemical potentials shows that the attraction is 
just a little bit to small for the occurence of the A's and one gets a pattern 
familiar from the RH-treatment (see Fig. 11). If one reproduces the same 
A-potential depths as in the RH-scheme, the A~ is still preferred at lower 
densities, but hyperons play now a more significant role due to their lower 
repulsion and dominate the high density region (see Fig. 12). With respect to 
the adjustment of the potential depths of the baryons, this RHF-model cor- 
responds exactly to the treatment within the RH-scheme. As in the case of 
nucleons, the different assumptions about the behaviour of the A-Dirac mass 
in matter influence the EOS only weakly. For instance, in the SU(6)-scheme 
for the hyperons the pressure increases only slightly by going from RHF 8 via 
RHF 9 to RHF 1 . More severe are changes of the repulsion of the baryons. 
For example, one obtains by increasing the hyperon repulsion by going from 
RHF 1 to RHF 4 {xi_jh = 0.8) pressure increases of ~150 MeV/fm^, compared 
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with ~ 30 MeV/fm^ in the case before with constant hyperon couphngs [|T3 
In conclusion we have constructed and discussed EOSs of NSM, which, 
in contrast to former investigations, are tight to the outcome of micro- 
scopic RBHF-calculations. The extrapolation to higher densities, where 
more baryons participate, was performed in the RH- and RHF-scheme. For 
the density domain of NSM, where the EOS is (completely) unknown, we 
were confronted as in all treatments with a complex composition with the 
necessity to invoke assumptions about the behaviour, especially with respect 
to the density dependence of the Dirac masses, controlled by an additional 
parameter. In this context we would like to emphasize again that so far 
almost all relativistic models have considered the rather special case of a 
RH-EOS based on a large Dirac mass at saturation (~ .79 rrifyf). With re- 
spect to the couplings of the hyperons and A's large uncertainties exist. For 
the hyperons the relative ratios of the cr- and u;- couplings were fixed by 
the potential depths in nuclear matter. (An absolute fixing from hypernuclei 
data is not possible at present For the A-coupling we invoked several 

assumptions, which take into account the smaller decrease of the A-mass in 
matter. 

The case of the RH-treatment resembles in its basic features to earlier 
investigations in this framework, which were based on more or less phe- 
nomenological treatments. Their distinguishing features are a high hyperon 
content (high strangeness), stiffer EOSs, and a strong A-suppression. In 
the RHF-scheme the A's play in general due to their smaller (yfp-coupling an 
important role, and the hyperons are more suppressed. The detailed com- 
position and stiffness of the EOS depends on the Dirac mass of the A, for 
which we selected four different choices. 

We are aware of the fact that the described uncertainties open the gates to 
a realm of options, some of which might be rather unfamiliar in comparison 
with standard treatments. Of course we could not explore and present all the 
hypotheses, but we have tried to select illustrative samples, from which the 
trends can be extracted (for more details, see Ref. [l^). The consequences 
for NSs will be discussed in the following sections. 
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Ill Neutron star properties 



a) General considerations 

The theoretical description of a neutron star is governed by the following 
conditions, which can be roughly estimated, for instance, on the basis of the 
Fermi gas or the hard core model for the EOS: i) General Relativity has to be 
taken into account for the determination of the gross properties of a star with 
approximately one solar mass and a radius R of approximately 10 km, 
since the relativistic effects (change of the metric etc.) for such objects are 
of the order 

ii) For the EOSs one can use relativistic treatments within a Minkowski 
metric since the spacing of baryons in the star {R = tqA^^^; tq ~ 0.5 fm, A ~ 
10^^) is of the order of 10~^^. iii) As explained before, due to the high 
densities one should favour relativistic EOSs for NSM. 

Roughly one can cast the EOSs in terms of stiff or soft equations. The 
"stiff" equations give a maximum mass (Oppenheimer-Volkoff limit Mov) of 
about twice the solar mass, and limiting rotational periods larger than 1.5 ms. 
For "soft" EOSs the estimates are Mqv = 1-5 Mq and Pmm = 0.6ms. Hence, 
the detection of sub-milliseond or heavy pulsars could discriminate the (the- 
oretical) EOSs. Furthermore the millisecond pulsars, explained according to 
the present understanding within the so-called recycling model, are very fast 
rotating objects (at birth at least 0.5 ms) Consequently it is important 
to include the (fast) rotation into the theoretical treatment of NSs. 

b) Theoretical treatment of rotating and deformed stars 

Due to these conditions one is faced with the problem to determine Ein- 
stein's curvature tensor G^u for a massive star {R^^.g^i,, and R denote the 
Ricci tensor, metric tensor, and Ricci scalar, respectively). 

G^, = R^, - ^g^^R = SvrT^, (e, P(e)) . (III.2) 

A necessary ingredient for solving (III. 2) is the energy-momentum tensor 
density T^y, for which knowledge of the (relativistic) EOS, i.e. pressure P as 
function of the energy density e is necessary. 
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For a spherically symmetric and static star, the metric has the famous 
Schwarzschild form {G = c = 1): 

ds^ = - e^'i'^'^dt' + e^^^'^dr^ + r\de^ + sin^ed<p') , (III.3) 

where the metric functions are given by: 

e2AM = (l-7(r))-\ (III.4) 

e2<^('^)=e-2^M = (l-7(r)) for r > R^tar , (III.5) 

with 

r ^ r<Rs 

7(r)= J, (III.6) 
y AM^ r>R, 

Einstein's equations for a static star reduce then to the familiar Tolman- 
Oppenheimer-Volkoff equation (TOV) ^, 0]: 

= - 1 (e(r) + P(r)) (M(r) + AiTr^P{r)) e-^^^'^ , (III.7) 

dr 

where the gravitational mass M(r) contained in a sphere with radius r is 
determined via the energy-density e(r) by: 



M(r) =A7T [ e{r)r^dr . (III.8) 
Jo 



The metric function 0(r) obeys the differential equation 

d(h 1 dP 



dr e(r) + P(r) ' 
with the boundary condition 



(III.9) 



<P{r = Rs) = ^Hl--f{R,)) . (III.IO) 

For a given EOS i.e. -P(e), one can now solve the TOV equation by integrating 
them for a given central energy density ec from the star's centre to the star's 
radius, defined by P{Rs) = 0. 
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More complicated is the case of rotating stars, where due to the rotation 
changes occur in the pressure, energy density, etc. The energy-momentum 
density tensor T^,^ takes the form {g^^u^Ut, = — 1) @, 0; EH' E31- 



T. 



+ AT, 



(III.ll) 



with 



rpO 



{e + P)u^u, + Pg^, , (III.12) 

AT^, = (Ae + AP)u^u, + APg^, . (III. 13) 

P, e, and p are quantities in a local inertial frame comoving with the fluid 
at the instant of measurement. For the rotationally deformed, axially- 
symmetric configurations one assumes a multipole expansion up to second 
order (P2 denotes the Legrendre polynomial): 



AP 
Ae 

Ap 



{e + P){po + P2P2{cos9)) 
dP ' 



AP 



dp 
dP 



(III.14) 
(III.15) 

(III.16) 



For the rotating and deformed star with the rotational frequency Q one has 
now to deal with a generalized Schwarzschild metric, given by 



(III.17) 



Here, uj{r) denotes the angular velocity of the local inertial frame, which - 
due to the dragging of the local system - is proportional to Q. 

The metric functions of Eq. (|III.17D which correspond to stationary rota- 
tion and axial symmetry with respect to the axis of rotation are expanded 
up to second order as (independent of and t): 



^2u{r,e,n) 
^2iiir,e,n) 



e^<t>{r) [1 + 2 (/io(r, n) + /i2(r, fi)P2(cos0))] 
sin^ ^ [1 + 2 (f2(r, ^]) - /i2(r, fi)) P2(cos t 



[1 + 2 {V2{r, Q) - /i2(r, Q)) P2(cos 

p2A(r) 



^ ^2 mo(r, il)G + m2(r, n)P2(cos( 
r 1 — 7(r) 



(III.18) 
(III.19) 
(III.20) 

(III.21) 
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The angular velocity in the local inertial frame is determined by the differ- 
ential equation 

A. ( r^a{r)^] + Ar'^^u;{r) = , r < , (III.22) 
dr \ dr J dr 

where uj{r) is regular for r = with ^ = 0. a(r) abbreviates 

a(r) = e~'^("Vl -7(0 • (III.23) 
Outside the star u{r, Q) is given by: 

cj(r, n) = n- ^J{n) , r> r, . (111.24) 

The total angular momentum is defined by: 

From the last two equations one obtains then an angular frequency f2 as a 
function of central angular velocity Uc = uj{r = 0) (starting value for the 
iteration) : 

Q{uj,) = uj{Rs) + J(^^) . (III.26) 

Due to the linearity of Eq. ( |111.22| ) for u){r) new values for uj{r) emerge simply 
by rescaling of Uc- The momentum of inertia, defined by / = is given by 
(a(i?.) = 1): 



/ =: ——^ = — / dvr* , =— — e"'^ . III.27 
n 3 Jo y/1 - 7(r) n ^ ^ 

Relativistic changes from the Newtonian value are caused by the dragging 
of the local systems, i.e. u/Q, the redshift (e"*^), and the space-curvature 
^(1 — 7(^))~^''^y For slowly rotating stars with low masses, one can neglect 
the dragging 1) and rotational deformations, but we would like to em- 

phasize that the described treatment is not restricted by low masses and/or 
slow rotations. 
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If one has determined ui{r), one solves in the next step the coupled mass 
monopole equations {i = 0) for mo,po ( = monopole pressure perturbation) 
and ho (for details, see Refs. 0). The quadrupole distortions h2 and V2 
{i = 2) determine the star's shape (see Refs. After the determination 

of the distortion functions, one can express the surfaces of constant density, 



the star's eccentricity ||50|, and the mass quadrupole moment of the star as 
follows: 

r{e) = r + Ur) + {Ur)+r[v2{r)-h2{r)]}P2{cose) , (III.28) 



A2 denotes an integration constant and is defined by {i = 0, 2): 



e. = -P.(e + p)(^) . (111.31) 



c) Stability criteria 

An intriguing problem in the physics of NSs is the question, whether the 
used EOSs are in accordance with the observed data, so supplying a test for 
the theoretical and partly speculative EOSs for highly compressed matter. 
Unfortunately, as far as gross properties (radii, masses) of NSs are concerned 
most of the nonrelativistic and relativistic EOSs are able to reproduce these 
gross properties (For a test of 25 different EOSs, see Refs. ^ ^ 0]). 

A more decisive criterion may be the stability of a star against rotation. 
Since no trivial stability criteria are known for rotating configurations in 
general relativity, we consider first Kepler's criterion, which sets an absolute 
upper limit on a star's rotation. The resulting Kepler frequency, fix, beyond 
which instability sets in due to mass shedding at the equator, is given for the 
generalized Schwarzschild metric as solution of {tp' = etc.) @, |^, ^ 
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with 



+ 



(III.33) 

to be evaluated at the equator. V denotes the orbital velocity of a comoving 
observer at the equator relative to a locally non-rotating observer. Neglection 
of the distortions (/i2, f 2 << 1) and of the dragging of the local inertial frames 
{iv = Q) gives 



eq 




leq 



Req^c 



or 



leq 



with 




(III.34) 



^111.35) 



The Newtonian result is recovered by using a flat space-time-geometry. The 
Kepler frequency of the heaviest neutron star can be obtained from the mass 
and radius of the most massive nonrotating neutron star as 



K 




(III.36) 



which has the advantage that one needs only the input from a static non- 
rotating star model. 

As mentioned before the Kepler criterion gives only an upper limit. Gravitational- 
wave reaction instabilities of a rotating star are likely to lower the star's 
maximum rotational frequency below VLk- Since the theory is rather lengthy, 
we refer for details to Refs. ^ ^ |5l], The critical frequency for a 
particular instability mode (m = 2, 3 . . .) is given by: 



T, 



g,m 



Tr 



(III.37) 



where u denotes the shear viscosity, depending on temperature T. The ex- 
pressions for the damping time scales r for gravitational radiation reactions 
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(Tg,m,), viscous damping (T^^m), and the surface mode Um can be found in 
Refs.[g0,|§0,g. 

d) Further quantities 

For completeness we also give the expressions for the redshifts, the injec- 
tion energy, and the stability parameter. 

The frequency shifts of light emitted at the equator in backward (b) and 
forward direction (f ) is given by 0, ^ : 

For the redshift at the pole one gets: 

Zp{n) = e-"^""^ - 1 . (111.39) 
The so-called injection energy is defined as: 

p^^^^^2umi = I . (111.40) 

For discussing the stability of rotating stars it is useful to define the stability 
parameter: 

m) := ^ . (III.41) 

where T denotes the rotational and W the gravitational energy of the star. 

e) Results (gross properties) 

As a first example for the influence of the rotation we show in Figs. 13 and 
14 the NS-mass versus central energy density and the radius-mass relation 
for a composition with p,n,e~, and /i~ only, for nonrotating stars and stars 
rotating with their Kepler frequency. As expected and discussed before one 
reaches, due to the large stiffness of the EOS, the maximum mass at lower 
central densities. The rotation can increase the mass by more than half a 
solar mass. In this context one may remark that relativistic EOS of nuclear 
matter are in general stiffer than their nonrelativistic counterparts [0, 
If one includes now hyperons in the relativistic NSM-EOSs one obtains a 
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softening of the EOS, so that the gross properties of NSs may not differ 
too much from NS-calculations with nonrelativistic treatments with a pure 
nucleonic/leptonic composition. Only the protons concentration are not the 
same in both cases, since the different behaviour of the symmetry energy 
lowers the proton contribution in the nonrelativistic case, so suppressing the 
so-called direct Urca-process ||53|, |55 |. 

Next we turn to the more interesting case of NSM including more baryons. 
For EOSs in the relativistic Hartree scheme one obtains - as expected - still 
sufficiently large enough NS-masses, which decrease - as discussed before - 
for weaker hyperon couplings. This behaviour is exhibited in Figs. 15 and 
16. So far the results comply with the familiar pattern. A limit for the rel- 
ative hyperon couplings, in this scheme, is given by the SU(6) choice, since 
otherwise the NS-masses become too small. This result compares with our 
NS-calculation with the NSM-EOS TMl, where we obtained a maximum 
star mass of 1.56 M© (see Section II). Furthermore we find the expected pe- 
culiar behaviour for the universal coupling in the lower energy-density region 
(see discussion in Section II). This can also be seen in Table V, where the NS- 
properties for a fixed mass of 1.4 Mq are given. They show with exception 
of the universal coupling the expected behaviour with increasing hyperon- 
couplings. For comparison we include also the outcome for nonrelativistic 
calculations without hyperons [^], which demonstrates clearly the softening 
of the relativistic EOS due to the hyperons (see also Figs. 13,15,16), which 
become even softer than the nonrelativistic EOS without hyperons. 

More interesting are the RHF-EOSs. In the case of using the hyperon 
couplings of the RH-treatment and Qn = one obtains for rather low 
densities an onset of A~ and hyperons, where the A~ and the S~ play a 
decisive role (neglecting the A~ does not change the EOS significantly, since 
then the S~ occur earlier |l^), and the EOS becomes rather soft for lower 
densities, but for higher densities the pressure rises again more strongly. As a 
consequence the gravitational NS-mass as a function of central energy density 
should show a flat plateau before it rises again in the standard pattern, but 
at higher densities than in the RH-EOS case. Furthermore one expects 
relatively low star masses due to the soft EOS. 

This behaviour is exhibited in Fig. 17, where the gravitational NS-mass 
is given as function of the central density. Reasonable star masses demand 
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larger hyperon couplings. The flat behaviour in the lower domain resembles 
- as one could have expected - to the case of rather soft EOSs with low 



incompressibility (see, for instance, the BCK-EOS [42], where one obtains 



maximal NS-masses of approcimately 1 Mq In Fig. 18 we exhibit the 

limiting Kepler frequencies. As expected the softer RHF-EOSs permit lower 
rotational periods. However inclusion of the gravitational-radiation instabil- 
ities shows that the Kepler frequency is only an upper limit on the critical 
frequency and the limiting periods increase of about 30 % @, ^ |10| (for more 



details, see Ref. 0]). If one uses RHF-EOSs, where the ratio of hyperon cou- 
plings is adjusted in the RHF-scheme (see Section II), the hyperons are not so 
easily produced and the resulting EOSs become stiffer (see Section II). In the 
presentation we will restrict ourselves to the interesting cases of the smallest 
hyperon couplings compatible with star masses around 1.5Mq. The resulting 
star masses as function of the central energy density are shown in Figs. 19 and 
20. By comparison with Fig. 17 we infer that these improved EOSs, where 
the hyperon potential depths are treated correctly, lead to masses, which are 
better in accordance with the data, since increase of the hyperon couplings 
and rotation lead to even higher star masses. EOSs with weak A-repulsion 
are rather soft and give only for strong hyperon couplings and high rotation 
frequencies mass values around 1.5 Mq (see Figs. 19,20). 

One may illuminate the situation further by comparing the properties for 
a typical NS of 1.4 Mq (see Tabs. V,VI). By the given arguments one should 
obtain increasing central densities by going from the stiff RH-EOS (RHl) 
to the softest RHF-EOSs with low Dirac masses for both the deltas and 
hyperons (RHF12,13). For the properly adjusted RHF-EOSs (RHF8,9) the 
star parameters differ not significantly. Here the higher A-masses suppress 
the influence of the A's below e ~ 500 MeV/fm^ and one obtains sufficient 
pressure to reach a mass of 1.4 Mq earlier (see also Fig. 19). Decrease of 
the A-mass according to universal coupling of the A's gives an increase of 
the central energy by a factor 2 (RHFl compared with RHl). Finally we 
illustrate in Table VII the influence of rotation. Shown are the results for 
a NS with the same baryon number. The mass changes due to rotation are 
relatively small for constant baryon number (for fixed central energy density 
see Refs [4,7,24]). 

f) Cooling properties 
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Another decisive test for the EOS may be the coohng history of a NS. In 
Ref. 1^ we have aheady calculated and discussed the thermal evolution of 



various models with different EOSs and different involved processes. Here we 
show and compare the cooling tracks of NS models constructed for the three 
EOSs RHl, RHFl, and RHF8, as well as the RHF-EOS without hyperons 
(see Fig. 22). During about the first million years the NS cools mainly 
by emission of neutrinos. One classifies the neutrino processes into slow and 
enhanced ones, according to whether two or only one baryon is participating. 
Enhanced processes cause a temperature inversion in young stars, i.e. the 
interior of the star becomes much cooler than the crust. Depending on the 
crust thickness the cooling wave formed by the temperature gradient reaches 
the surface and causes the sharp decrease of the surface temperature after 
about 30 years (see the three broken curves in Fig. 22). 

For the EOSs considered here the only possible enhanced cooling processes 



are the nucleon [36, 37 



n^p + l- + ui (III.42) 
and the hyperon direct Urea processes 



^- ^A + r + iJi (III.43) 

and 

A ^ p + /- + z7z , (III.44) 

as well as their inverse reactions. Here, I denotes electrons and muons. The 
nucleon direct Urea process is only possible, if the proton fraction exceeds 
some critical value of about 11 % for a pure nucleonic/electron composition of 
the neutron star matter, since otherwise energy and momentum conservation 
cannot be fulfilled simultanously. If hyperons and muons are taken into 
account this value rises slightly to approximately 13 %. Similar constraints 
have to be considered for the hyperon direct Urea processes. It is obvious 
that the resulting thermal evolution depends strongly on the EOS. It seems 
to be a general feature that non-relativistic EOSs have protron fractions 



below this critical value [^, whereas relativistic EOSs allow for the nucleon 
direct Urea process. The critical masses above which the hyperon direct Urea 
processes are possible are approximately equal to 1.3 Mq for all three EOSs 
studied in this section. This seems to be surprising, since A and appear 
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beyond n ~ 0.7 fm~'^ in the case of RHFl (see Fig. 9), and already beyond 
n ~ 0.3 fm~^ in the cases of RHl and RHF8 (see Figs. 6 and 11). This higher 
threshold density is however compensated by the smaller incompressibility 
of RHFl compared to RHl and RHF8 (see Sect. II). The used slow neutrino 
processes, as well as the processes in the crust of the NS, are discussed in 
greater detail in Ref. |56|. 



The cooling behaviour of a NS is also influenced by the appearance of 
superfiuid phases. If neutrons or protons become superfiuid the neutrino 
emissivity of the nucleonic processes, the thermal conductivity and the heat 
capacity are reduced by an approximately exponential factor exp(— A/ZcT), 
where A denotes the gap energy (see Table IV of Ref. for the used gap 
energies) . 

The observational data are described in Ref. |^8[ (see Table 2 in Ref. ||58|| ). 



The obtained effective surface temperature depends crucially on whether a 
magnetized hydrogen atmosphere is used or not. Since the photon flux, 
measured solely in the X-ray energy band, does not allow to determine what 
atmosphere one should use, we consider both the blackbody model (solid 
error bars in Fig. 22) and the hydrogen-atmosphere model (dashed error 
bars). The plotted errors represent the 3a error range due to the small 
photon fluxes. 

All models exhibit enhanced cooling via the nucleon direct Urea process. 
However this process is suppressed below the critical temperature for the 
superfiuid phase transition. Since the process ( [111.43| ) is not suppressed by 
superfluidity, the surface temperature of these models (see broken curves) 
is much smaller than the one of the model without hyperons (solid curve). 
The observed data can almost perfectly be described by the latter model, 
provided one assumes that the pulsars have no hydrogen atmosphere (except 
PSR 1055-52, which could be explained by internal heating; see, e.g. Refs. 

|60|). However, if some of the pulsars prove to have a hydrogen atmo- 
sphere, these models seem to be too hot. Whether the observed pulsars have 
a hydrogen atmosphere could be decided if one considers multiwavelength 
observations, as suggested by Pavlov et al. [0. With respect to cooling 



properties one seems to get along with simpler relativistic EOSs without hy- 
perons. However the weaker temperature drop in this case is caused by the 
superfluidity of nucleons, which cannot be included, at present, for the other 
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baryons. Inclusion of this effect for hyperons may shift the curves towards 
the observed values, since the process ( [ill.43| ) would be suppressed, too. 

Please note that we have considered only some of the possibilities of 
neutron star cooling. Additional processes as internal heating [ pU| , |BD[ or 
intermediate neutrino processes may yield different cooling tracks. This 



is also true for the effect of accreted atmospheres investigated in Refs. ^3 
61. 



IV Summary 

The goal of this investigation is to incorporate "parameter-free" microscopic 
relativistic Brueckner-Hartree-Fock calculations of nuclear matter in the in- 
vestigation of neutron star matter. In a first step we extended the RBHF- 
theory of asymmetric nuclear matter to neutron star matter, consisting of 
neutrons, protons and leptons, which has to obey the constraints of charge 
neutrality and generalized /3-equilibrium. Since for higher densities more 
baryons (hyperons etc.) have to be included, for which, at present, mi- 
croscopic RBHF-calculations are not feasible, we extended the scheme by 
utilizing either the relativistic Hartree- or Hartree-Fock-approximation, in 
which the other baryons can be incorporated. The coupling constants of these 
schemes were adjusted in the nucleonic sector to the outcome of the RBHF- 
calculations of NSM around saturation densities. In this manner we obtained 
a good description of NSM near saturation, which is essential for lighter neu- 
tron stars and also important for heavier stars. As long as one restricts the 
composition to p,n,e~, and /i" only, one obtains in this framework, due to 
the rather stiff EOS, Oppenheimer-Volkoff star masses around 2.2 M© and 
minimum rotation frequencies slightly above 1 ms. However if one incorpo- 
rates more baryons in the scheme one obtains a considerable softening of the 
NSM-EOS, which leads inevitably, as in the case of realistic phenomenologi- 
cal parametrizations of the nuclear Lagrangian, to negative Dirac masses for 
the nucleons in NSM. This drawback is the result of the necessity to repro- 
duce a rather low nucleon Dirac mass at saturation, which leads to a peculiar 
feature of the mean field approximation, namely that such Lagrangians cause 
a strong decrease of the effective cr-mass at higher densities, resulting in a 
steep decrease of the Dirac masses. Since the behaviour of the Dirac masses 
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at high densities arc (completely) unknown in NSM, we extrapolated the 
Lagrangian in these domains by a slightly changed dynamics, where the de- 
crease of the Dirac mass is not so severe, which is also in accordance with 
RBHF-results in nuclear matter. The resulting pressure changes are rather 
small. Within the RH-framework we obtained then compositions of NSM 
which are more or less in accordance with former investigations, which use a 
priori in the whole domain Lagrangians with increasing effective a-masses. 
For the hyperons we used different coupling strengths but the ratios of the 
a — a;-couplings were fixed by utilizing the potential depths of the hyperons 
in nuclear matter. Due to large p-coupling in the RH-approximation the A's 
are negligible in this framework. The Oppenheimer-Volkoff (OV) star mass 
reaches from 1.5 Mq - 2.4 depending on the hyperon couplings and the 
rotation frequencies. More complicated is the situation for RHF-EOSs. Here 
one is confronted with smaller p-couplings, which favours A's in the com- 
position, and the fact that the potential depths of the hyperons in nuclear 
matter are solely determined by their Hartree contributions. As discussed in 
details in Section II, the resulting compositions depend now strongly on both 
the assumptions about the hyperon couplings and the assumed behaviour of 
the A-Dirac masses in NSM. As long as one assumes for the A's the same 
coupling as for the nucleons, the A's play a dominant role in the composi- 
tion. The EOSs are relatively soft and the minimal Oppenheimer-Volkoff 
mass is around 1.5 M0 for weak hyperon couplings. This mass drops even to 
1.3 Mq for strong hyperon couplings if one decreases the relative A~-strength 
generally to 0.625 . Increase of the A-repulsion gives compositions not so 
different from the RH-approximations with minimum OV-masses of 1.6 M©. 

In general we can conclude in accordance with earlier findings that the 
EOSs based on hadronic theories of matter are capable of accomodating the 
gross properties as well as rotational periods of all pulsars known to date. 
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Table captions 



Table I: Saturation properties of infinite nuclear matter. RBHF-A and 
RBHF-i? denote the RBHF-results for the Brockmann potentials A 
and B, calculated in the full Dirac space (for details, see Ref. For 



comparison we added some results of more phenomenological relativis- 
tic Hartree calculations (NLl, NL-SH) ||23| and phenomenological non- 
relativistic calculations. SkM* and S 111 denote two well known Skyrme 
forces, TF96 is a recent Thomas-Fermi calculation [24]. Also included 



are two nonrelativistic microscopic variational calculations ||25 



Table II: Parametrizations of the RH- and RHF-Lagrangian adjusted to the 
RBHF-calculations. For the masses the following values were selected 
(MeV): TTiN = 939, m„ = 550, = 738, = 138, = 770 
(^f^ = 1.00265 f^/4:iT = 0.08). The parametrizations are labelled as 
follows: RHA = RBHA etc. for the potential A (effective mass at the 
Fermi surface m = 617.8 MeV (A); 621.8 MeV (B)). 

Table III: EOSs for the different density regions of a NS (1 MeV/fm^ corre- 
sponds to 1.783 xlO^^g/cm^). For the density region above 20 MeV/fm^ 
we use the parametrizations of the RBHF-calculations in the frame of 
the RH- and RHF-approximation, respectively (OBE-potential B of 
Brockmann and Machleidt). The parameters for the nucleonic sector 
are given in Table II. For the population of the more massive baryons 
the 12 lowest lying ones are allowed for. 

Table IV: Relative coupling strengths of the hyperons in the different ap- 
proximations (see text). The ratios x^h/xuih are adjusted to the bind- 
ing energy of the hyperon in nuclear matter. For RHFll - RHF14 the 
RH-couplings are used. The universal coupling is defined hj = gn = 
g^. The EOSs for RHF7, RHFll and RHF12 are too soft for obtain- 
ing OV-masses around 1.5 Mq (see Fig. 17). In these models larger 
hyperon couplings are needed, for instance, by going from RHF7 to 
RHFIO Mov increases approximately to 1.5 Mq (at Kepler frequency). 

Table V: Comparison of the properties of a static, spherical NS of mass 
1.4 M0 for different RH-models. edPc) denotes the central energy 
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density (pressure). The amu mass Ma minus the gravitational mass 
Mq is effectively the binding energy liberated when the NS is formed. 
R denotes the star's radius, Ac stands for the stellar crust using 2.4 
X 10^"^ g cm~'^ as the boundary, I denotes the moment of inertia and z 
the surface redshift. For a comparison we included two nonrelativistic 
models without hyperons (see Table I) |p3| . 



Table VI: Comparison of the properties of a static, spherical NS of mass 
M = IAMq for different RHF models. Labels as in Table V. 

Table VII: Properties of rotating neutron star models of rotational period 
P = 1.4 ms and the same baryon number as the nonrotating star with 
M = IAMq, calculated for different EOSs. The entries are: central 
energy density ec, equatorial and polar radii, -Req and -Rp, respectively; 
moment of inertia /; stability parameter, t; injection energy f3; red- 
shift of the pole, z^, eccentricity, e quadrupole moment, 11. The 



gravitational mass increase due to rotation is rather small. 
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Figure captions 



Fig. 1: EOSs for asymmetric nuclear matter. Compared are the EOSs for 
the Brockmann-Machleidt potential B for a fixed asymmetry 6 in the 
RBHF-approximation (RBHF-B: full curves) with the treatment in the 
RH (RHB: dotted curves) and RHF (RHFB 1: dashed curves) -scheme, 
respectively. 

Fig. 2: Comparison of the EOSs for neutron star matter composed of n, p, e~, 
and fi- (RBHF-B, RHB, RHFB 1; Brockmann potential B). The upper 
branches correspond to the case where myons are neglected. 

Fig. 3: Comparison of the nucleon/lepton composition of neutron star mat- 
ter (potential B). The upper branches correspond to the case without 
myons. 

Fig. 4: Comparison of the Dirac masses with zero charge: The full curves 
show the behaviour in the original scheme. The dotted curves give 
the extrapolation according to Eqs.( [il.lD| , [il.llD . The selected case is 
the RHF-approach with SU(6) couplings (RHFll, a = 0.02) (see Ta- 
ble IV). The other cases are rather similiar. The masses differ hardly 
for the different isospin states. 

Fig. 5: Relative baryon/lepton populations in the relativistic Hartree scheme 
for universal coupling of the hyperons and deltas (RHl; a = 0.015) 

Fig. 6: Relative baryon/lepton populations in the relativistic Hartree scheme 
for hyperon-SU(6)-coupling of the vector mesons (RHl; a = 0.015). 
The relative a-couplings of the hyperons are adjusted to the corre- 
sponding potential depths in nuclear matter. 

Fig. 7: Dependence of the EOS on the choice of the relative a-meson- 
hyperon coupling Xo-h (relativistic Hartree EOS, potential B, a = 
0.015). Compared are the choices x^h = 0.7,0.8 with the so-called 
SU(6) parametrization and the universal coupling (see Table IV). For 
the first two cases the (larger) hyperon-vector couphngs are adjusted 
to the potential depths in infinite matter (for x^^h = 0.9,^^^^^ > 1). 
For the SU(6) parametrization the vector meson couplings are fixed by 
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the quark picture and the (smaller) a-couplings are adjusted to the 
potential depths {xo-h < x^h)- Due to the stronger a-coupling in the 
universal case one gets for lower densities a softer EOS, otherwise one 
confirms the expectation that smaller couplings softens the EOS (see 
text). 

Fig. 8: Relative baryon/lepton populations in the RHF-approximation for 
relative hyperon couplings taken from the RH-approximation in the 
SU(6) case (RHFll; a = 0.02) 

Fig. 9: Relative baryon/lepton population in the RHF-approximation (RHFl; 
a = 0.02). The cr-couplings are adjusted to the hyperon potential 
depths (see Table IV). Increase of the hyperon coupling would give a 
later onset of the hyperons, for instance, A occurs at p ~ 0.8 fm~^. 

Fig. 10: Relative baryon/lepton population for NSM in the RHF-approximation 
with reduced A-coupling {x^h = SU(6), x^h adjusted, (RHF7; a = 
0.02), xaa = xan = 0.625 [0). The onset of the As/hyperons is now 



shifted to higher/lower densities than for = gN- 

Fig. 11: Relative baryon/lepton population for NSM in the RHF-approximation 
(RHF8) with fixed x^-a = 0.625 and x^ja = 1- The hyperon couplings 
are adjusted to the hyperon-potential depths. 

Fig. 12: Relative baryon/lepton population for NSM in the RHF-approximation 
(RHF9) with adjusted delta- and hyperon couplings to the potential 
depths. 

Fig. 13: Gravitational star mass (in units of solar Mass Mq) as a function 
of central energy density for star models constructed from EOSs with 
p,n,e~, and The upper curve corresponds to (deformed) stars 
rotating at their Kepler frequency. The RHF- and RH-curves are very 
close, since the pressure differs not much. 

Fig. 14: Neutron star radius versus mass. Shown are sequences of stars 
rotating at their Kepler frequencies and at zero frequency {p,n,e~, 
and fi~ only). The minimal periods are ~ 1 ms for M = 1.5M0. 
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Fig. 15: Dependence of the gravitational star mass on the relative meson- 
hyperon coupling. The non rotating, spherical star families are given 
as function of the central energy density (RH-EOSs, a = 0.015; all 
baryons included). The maximum star mass increases with stronger 
hyperon couplings. The peculiar behaviour for the universal coupling 
for lower densities is explainable by the stronger attraction (see Figs. 5,7 
and text). 

Fig. 16: Increase of the gravitational star mass for rotating stars: Shown 
are the star families as in Fig. 15, but now for deformed stars rotating 
at their Kepler frequency. The radii are approximately 12 (15) km for 
nonrotating (rotating) stars at Mj^s = 1-5 Mq. 

Fig. 17: Gravitational NS mass for nonrotating stars versus central energy 
density for RHF-EOSs with universal couplings and couplings from the 
RH-treatment. The peculiar behaviour for low densities corresponds 
to the early onset of baryons in this model (see text). At their Kepler 
frequency the OV-mass ranges from ~ 1.95 Mq (universal) to ~ 1.5 Mq 
(SU(6)-model); for M = I.IMq the radius stretches in the range 9 km 
(nonrotating) to 11 km (Kepler rotating). 

Fig. 18: Limiting rotational Kepler periods of pulsars versus NS mass for 
the models described in Figs. 16 and 17. The shaded area covers the 
range of observed periods and masses . 



Fig. 19: Gravitational NS mass for spherical nonrotating stars versus cen- 
tral energy density. Compared are the star families for EOSs in the 
RHF-scheme for different A-couplings for the weakest hyperon cou- 
pling compatible with Mqv ~ 1-5 Mq. The ratio of the hyperon cou- 
plings is adjusted to the hyperon potential depths in nuclear matter. 

Fig. 20: Gravitational NS mass for nonspherical stars rotating at their Ke- 
pler frequency versus central energy density. Compared are the same 
EOSs as in Fig. 19. 

Fig. 21: Limiting rotational Kepler periods of pulsars versus NS mass for 
the models described in Figs. 19, 20. Inclusion of gravitation-radiation 
instabilities increases the limiting period by approximately 30% (for 
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M = 1.5 Mq) |18|. The shaded area covers the range of observed 
periods and masses |EB . 



Fig. 22: Coohng of M = IAMq models for different EOSs. The surface tem- 
perature obtained with a blackbody- (magnetic hydrogen-) atmosphere 
are marked with sohd (dashed) error bars labeld by the respective pul- 
sar's position. 
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TABLE I 





E/A 


Poo 




J 




(MeV) 


(fm-3) 


(MeV) 


(MeV) 


RBHF-A 


-16.49 


0.174 


280 


34.4 


RBHF-B 


-15.73 


0.172 


249 


32.8 


NLl 


-16.4 


0.152 


212 


43.5 


NL-SH 


-16.3 


0.146 


356 


36.1 


SkM* 


-15.8 


0.160 


216 


30.0 


SIII 


-15.9 


0.145 


355 


28.2 


TF96 


-16.24 


0.161 


234 


33 


WUU 


-15.5 


0.175 


202 


30 


WUT 


-16.6 


0.157 


261 
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TABLE II 





9a 




9p 


10^ X hN 


10^ X Cjv 


fp/9p 


RHA 


9.58096 


10.67698 


3.81003 


3.333665 


-3.52365 




RHFAl 


9.28353 


8.37378 


2.10082 


3.333689 


-2.15239 




RHFA2 


9.24268 


8.25548 


2.19809 


2.96514 


-2.68614 


3.7 


RHFA3 


9.16665 


8.07540 


2.37987 


1.95524 


-2.36335 


6.6 


RHB 


9.59169 


10.68084 


3.66541 


3.62616 


-4.17140 




RHFBl 


9.36839 


8.40466 


1.77326 


3.74354 


-3.18456 




RHFB2 


9.33266 


8.32154 


1.86078 


3.44306 


-3.46261 


3.7 


RHFB3 


9.26782 


8.19391 


2.02216 


2.67292 


-3.18198 


6.6 
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TABLE III 



EOS 


Energy-density range 


composition 




(MeV fm^) 




HW^ 


e < 0.6 


Crystalline; light metals, electron gas 


NV'' 


0.6 < e < 20 


Metals, relativistic electron gas 


RHB 


20 < e 


p,n,e-,/i- 


RHFBl 


20 < e 


p,n,e-,iJ,- 


RHB 


20 < e 


p, n, A, A's, e-,fi- 


RHFBl 


20 < e 


p,n,A,S±'°S°'-, A's,e-,/i- 



taken from Ref. [IS 



taken from Ref. 0] 
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TABLE IV 



Approximation 














RH1(SU(6)) 


0.676 


2/3 


0.342 


1/3 


1 


1 


RH2 


0.7 


0.77 


0.7 


0.78 


1 


1 


RH3 


0.8 


0.9 


0.8 


0.91 


1 


1 


RH4 


0.9 


1.03 


0.9 


1.04 


1 


1 


RH5(U) 


1 


1 


1 


1 


1 


1 


RHF1(SU(6)) 


0.44 


2/3 


0.26 


1/3 


1 


1 


RHF2 


0.36 


0.5 


0.35 


0.5 


1 


1 


RHF3 


0.46 


0.7 


0.45 


0.7 


1 


1 


RHF4 


0.51 


0.8 


0.5 


0.8 


1 


1 


RHF5 


0.56 


0.9 


0.55 


0.9 


1 


1 


RHF6(U) 


1 


1 


1 


1 


1 


1 


RHF7(SU(6)) 


0.44 


2/3 


0.26 


1/3 


0.625 


0.625 


RHF8(SU(6)) 


0.44 


2/3 


0.26 


1/3 


0.625 


1 


RHF9(SU(6)) 


0.44 


1/3 


0.26 


1/3 


0.75 


1 


RHFIO 


0.51 


0.8 


0.5 


0.8 


0.625 


0.625 


RHF11(SU(6)) 


0.676 


2/3 


0.342 


1/3 


1 


1 


RHF12 


0.7 


0.77 


0.7 


0.78 


1 


1 


RHF13 


0.8 


0.9 


0.8 


0.91 


1 


1 


RHF14 


0.9 


1.03 


0.9 


1.04 


1 


1 
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TABLE V 



Quantity 


RH1(SU(6)) 


RH2 


RH3 


RH5(U) 


WUT 


TF 


Ce (10^^g/cm=^) 


8.1 


7.7 


7.2 


8.9 


12.12 


10.30 


Pc (lO-'^Myn/cm^) 


8.98 


8.62 


8.22 


11.79 


18.58 


14.54 


Mg/Mq 


1.4 


1.4 


1.4 


1.4 


1.4 


1.4 


(M^ - Mg)/Mo 


0.166 


0.166 


0.166 


0.163 


0.169 


0.173 


R(km) 


12.79 


12.84 


12.89 


12.21 


10.86 


11.37 


Ae (km) 


1.88 


1.89 


1.92 


1.67 


0.88 


1.08 


I (10^^ g cm2) 


16.14 


16.28 


16.46 


14.58 


12.30 


13.34 


z 


0.215 


0.214 


0.213 


0.229 


0.271 


0.254 



TABLE VI 



Quantity 


RHFl 


RHF6 


RHF8 


RHF9 


RHF12 


RHF13 


ec (lO^Vcm-'^) 


16.3 


16.8 


9.3 


10.1 


21.8 


17.8 


Pc (10^Myn/cm2) 


29.52 


33.59 


10.09 


11.91 


45.2 


35.47 


Mg/Mq 


1.4 


1.4 


1.4 


1.4 


1.4 


1.4 


{Ma - Mg)/Mq 


0.21 


0.25 


0.16 


0.16 


0.216 


0.219 


R(km) 


9.86 


9.41 


13.15 


12.8 


9.31 


9.59 


A, (km) 


0.8 


0.81 


2.4 


2.25 


0.88 


0.98 


I (10^^ g cm2) 


10.43 


9.86 


16.11 


15.16 


9.25 


9.80 


z 


0.312 


0.335 


0.208 


0.215 


0.341 


0.325 
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TABLE VII 



Quantity 


RHl 


RH5 


RHFl 


RHF8 


(MeV/fm^)) 


403.12 


479.40 


869.83 


438.85 


-Req (km) 


13.72 


12.95 


10.14 


14.49 


i?p (km) 


12.11 


11.68 


9.65 


12.6 


log //(gcm^) 


45.17 


45.13 


45.00 


45.19 


t 


0.036 


0.03 


0.016 


0.038 




0.658 


0.646 


0.572 


0.672 




0.232 


0.244 


0.323 


0.22 


e 


0.47 


0.43 


0.31 


0.49 


n (km^) 


7.27 


5.14 


1.51 


8.32 
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